AD-A070  793  OHIO  STATE  UNIV  COLUMBUS  ELECTROSCIENCE  LAB  F/6  17/9 

RADAR  TAR6ET  IDENTIFICATION. (U) 

MAY  79  0 M MOFFATT » K A SHUBERT.  E M KENNAUGH  F19628-77-C-0125 

UNCLASSIFIED  ESL-784677-2  RADC-TR-79-118  NL 


HLt  lAjn  HBA070793 


0.  H.  Moffatt 
K.  A.  Shubert 
E.  M.  Kennaugh 


Ifcis  report  has  been  reviewed  by  the  MDC  Inforoation  Office  (01)  end 
le  releasable  to  the  National  Technical  Iafornatloo  Service  (NTI8) . At  NTIS 
it  will  be  releasable  to  the  general  public.  Including  foreign  natlone. 

RADC-TR-79-118  has  been  reviewed  and  is  approved  for  publication. 


'■''mK'-d 

4. "l 


jo m p.  boss 

Acting  Chief,  Plans  Office 


nailing  list,  or  if  the  addressee  Is  no  longei 
please  notify  RADC  («BC)  Hanscon  APB  MA  01731 
nalntalnlng  a current  nailing  list. 

Do  not  return  this  copy.  Retain  or  destroy. 


p IQVPP 


UNCLASSIFIED 

SECUR|TX_QASSIFlCATION  OF  THIS  PAGE  (Wh on  P«C  Enur'd) 

F Ov  REPORT  DOCUMENTATION  PAgT 


REP  TOT  NUMBER 


12.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
3.  RECIPIENT'S  CATALOG  NUMBER 


1 4.  TITL6  (md  Subtitle) 


£ |3ADAR  ^ARGET  IDENTIFICATION  # 


5.  TYPE  OF  REPORT  A PERIOD  COVERED 

Final  technical  J^epStt  7~) 
i 1 Mar  77—  36  Sen  78,  S 

"^PERFORMING  Oka.  REPORVNUMBER 


17.  AUTHOWQ) 


f/D  FD-  M./faoffatt. 
U V |K.  A./Shubert7 


fANT  NUMBERC.; 


F19628-77-Ct6i25 


V PKHH3RMINU  ORGANIZATION  NAME  AND  AODRESS 

The  Ohio  State  University  ElectroScience 
Laboratory/Dept  of  Electrical  Engineering 
Columbus  OH  43212 

II.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

Deputy  for  Electronic  Technology  (RADC/EEC) 
Hanscom  AFB  MA  01731 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  « WORK  UNIT  NUMBERS 


IK.  BEport  DATE 
1 May  79  / 

' U.  NUMBER  OF  PAGES 


4.  MONITORING  AGENCY  NAME  a ADDRESSCIf  dllleront  tram  Controlling  Ollico)  I 15.  SECURITY  CLASS,  (ol  II 


IIS.  DISTRIBUTION  STATEMENT  (ol  Ihlo  Report) 


IS*.  OECLASSIFIC  ATI  ON/ DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


I 17.  DISTRIBUTION  STATEMENT  (of  the  obo  tract  entered  In  Block  20,  It  different  front  Report) 


IB.  SUPPLEMENTARY  NOTES 

RADC  Project  Engineer:  Ortho  E.  Kerr  (EEC) 


It.  KEY  WORDS  (Continue  on  rovoroo  elde  II  neceeee 

Radar  Resonances 

Target  Swept 

Identification  Frequency 

Natural 


end  Identify  by  block  number) 


10.  ABSTRACT  (Continue  on  rover ee  aide  If  neceaaery  end  Identify  by  block  number) 

This  report  presents  a review  of  the  radar  target  identification  problem  with 
primary  emphasis  on  practical  solutions  as  opposed  to  abstract  inverse 
scattering  formalisms.  A time  domain  viewpoint  is  stressed  throughout  because 
of  its  compactness,  ease  of  visualization  and  diagnostic  potential  and  because 
It  Is  a unifying  tool  for  various  electromagnetic  approximations.  A small 
finite  set  of  excitstlon  Invariant  target  descriptors,  the  coefficients  of  a 
linear,  homogeneous  difference  equation,  or  equivalently  the  complex  natural  _ 


on  FO"“ 

w • JAN  7* 


m n\- 


UNCLASSIFIED 


AFCnNITY  Classification  nTTm  B.nr  n».. 





SECURITY  CLASSIFICATION  OF  THIS  FAGEfRhwi  Dim  Bnlmnd) 


^-resonances,  are  justified  as  a discrimination  tool.  A prediction-correlation 
method  for  using  these  coefficients  in  identification  is  demonstrated.  It  is 
shown  that  substructure  features  of  the  targets,  in  this  case  vertical  and 
horizontal  stabilizers  of  aircraft,  possibly  can  be  used  as  a discrimination 
feature  and  that  the  required  interrogating  radar  signal  spectrum  for  this 
case  is  significantly  higher  than  that  required  for  previous  discrimination 
schemes  using  excitation  invariant  parameters.  In  this  specific  case,  however, 
the  target  parameters  do  not  have  a simple  physical  Interpretation. 
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PREFACE 

One  need  only  be  mildly  optimistic  to  forecast  that  future  gen- 
erations of  radar  systems  will  routinely  provide  a great  deal  of  target 
Information  beyond  the  usual  range,  elevation,  azimuth  and  velocity. 

At  the  very  least  these  new  systems  will  provide  parametric  target 
data  for  diagnostics  and  the  other  extreme  will  be  fully  automated 
target  Identification.  The  key  to  this  new  sophistication  will  be 
the  exploitation  of  broadband  complex  scattering  data  spanning  several 
octaves.  The  ElectroSclence  Laboratory  of  The  Ohio  State  University's 
Electrical  Engineering  Department  has  pioneered  the  interpretation 
and  utilization  of  multiple  frequency  complex  scattering  data  for 
numerous  purposes  including  radar  target  identification.  The  extrap- 
olation of  these  results  to  new  generations  of  radar  systems  is  the 
logical  next  step.  In  the  Interim  preceding  these  developments  it 
Is  meaningful  to  explore  what  can  be  accomplished  if  the  basic  broad- 
band premise  is  deliberately  distorted  to  more  nearly  accommodate 
existing  radar  systems.  Much  of  the  subject  matter  of  this  report 
is  directed  to  this  question. 
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EVALUATION 


The  experimental  and  theoretical  work  done  on  this  contract 
shows  that  aircraft  of  similar  size  and/or  different  models  of  the 
same  aircraft  (examples  are  F104  and  F104N)  can  be  uniquely  identified 
by  processing  a radar  signal  which  contains  a number  of  frequencies. 

The  techniques  developed  use  the  resonances  of  the  aircraft  which 
are  invariant  with  respect  to  the  aspect  angle,  however  their  amplitude 
do  vary  with  aspect.  The  resonances  are  translated  into  mathematical 
pole-zero  representations  for  each  aircraft  and  stored  as  data  in  a 
library  to  compare  with  the  data  received  from  an  unknown  aircraft. 

This  work  is  in  support  of  TPO  R2B.  Although  the  frequencies  used 
in  this  effort  are  lower  than  the  present  radars  use,  the  technique 
may  possibly  be  used  at  higher  frequencies  to  detect  smaller  resonant 
parts  of  aircraft  that  are  unique  to  each  aircraft. 

(fjk  r/^ 
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INTRODUCTION 


The  subject  matter  of  this  report  is  the  discrimination  of  radar 
targets;  specifically  aircraft  targets  of  the  jet  fighter  class.  Disc- 
rimination* rather  than  identification  is  used  here  to  denote  that 
the  techniques  to  be  described  distinguish  targets  wi th i n a given 
class.  Targets  not  in  the  class  can  at  best  be  categorized  as  being 
outside  the  class  and  possibly  looking  most  like  a particular  target 
in  the  class.  It  is  not  intended  to  imply  from  the  above  that  tar- 
gets of  disportionate  size  as,  e.g.,  bombers  and  light  aircraft,  can- 
not be  discriminated.  This  in  fact  is  a simple  discrimination  task. 
Conversely,  the  discrimination  of  targets  of  nearly  equal  size  and 
similar  shape  is  a most  difficult  task.  It  is  this  difficult  task 
which  is  addressed  in  this  report. 

Active  radar  interrogation,  regardless  of  the  particular  radar 
system,  basically  involves  the  interaction  of  electromagnetic  waves 
and  material  objects.  From  a system  viewpoint  the  radar  cross  section 
of  a target  is  a statistical  quantity.  However,  to  achieve  an  insight 
into  the  Identification  problem  it  is  necessary  to  return  to  a fun- 
damental deterministic  model.  The  scattered  electromagnetic  field 
from  an  object  is  controlled  by  the  parameters  listed  in  Table  I. 

Only  the  parameters  at  the  left  in  Table  I are  at  the  degression 
of  the  system  designers.  To  achieve  a target  detection  and  identi- 
fication capability  therefore  requires  a systematic  processing  of 
selected  parametric  variations  available  to  the  designer.  To  select 
the  appropriate  parameters,  however,  requires  a fundamental  under- 
standing of  how  each  of  the  target  characteristics  at  the  right  in 
Table  I affects  the  scattered  field.  It  is  interesting  to  observe 
that  the  above  remarks  could  just  as  easily  be  applied  to  radar  cross 
section  control  if  one  interchanges  the  controllable  and  noncontrollable 
parameters. 

♦Hereafter  the  words  identification  and  discrimination  are  used  inter- 
changeably. 
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Table  I 


Frequency 

Target  Size 

Time 

Target  Shape 

Radar  Polarization 

Target  Composition 

Pulse  Shape  and  Length 

Target  Orientation 

Pulse  Repetition  Frequency 

Target  Velocity 

Source-Receiver  Separation 

Component  Target  Motion 

Scattering  Parameters 


Radar  practitioners  will  immediately  observe  that  the  last  three 
items  on  the  right  in  Table  I,  orientation,  velocity  and  component 
motion,  can  and  are  exploited  as  a basis  for  target  identification. 

For  example,  a combination  velocity  pattern  and  heading  usually  in- 
dicates a particular  mission  and  consequently  a particular  vehicle. 

Component  motion  such  as  the  rotor  blades  on  helicopters  place  a par- 
ticular modulation  on  the  radar  signal.  Other  mechanisms  yield  modu- 
lations of  the  signal  or  perturbations  of  the  doppler  spectrum.  We 
do  not  pursue  such  mechanisms  because  they  are  not  fundamental  to 
the  electromagnetic  interactions  of  the  basic  target.  That  is,  a 
modulation  induced  by  a small  antenna  on  a spherical  object  may  in- 
deed be  recognizable  but  the  relationship  to  the  basic  spherical  ob- 
ject is  secondary.  By  the  same  token,  a given  complex  target  may 
exhibit  particular  modulation  or  doppler  characteristics  but  the  re- 
lationship to  the  object  is  secondary.  Such  characteristics  would 
also  appear  to  be  easily  camouflaged  or  distorted.  Conversely,  the 
identification  parameters  pursued  on  this  contract  are  fundamental 
to  either  the  overall  physical  properties  of  the  radar  targets  or 
a projected  subset  of  them.  The  use  of  identification;  friend,  foe  f 

or  neutral  (IFFN)  techniques  with  transponders  have  also  been  ignored. 

It  is  unlikely  that  unfriendly  targets  will  be  cooperative  and  respond  . 

to  interrogation  overtures.  Even  friendly  targets  may  find  it  unwise 

2 


to  respond  because  of  tactical  considerations.  The  identification 
procedures  pursued  on  this  contract  do  not  require  the  cooperation 
of  the  target. 

The  first  three  parameters  at  the  left  in  Table  I are  the  most 
fundamental  quantities  in  terms  of  eliciting  physical  properties  of 
the  target.  This  insight  has  been  achieved  via  the  impulse  response 
concept  * as  applied  to  electromagnetic  scattering  problems.  One 
of  the  first  things  that  became  apparent  was  that  the  gross  physical 
properties  (size,  shape  and  composition)  of  the  target  were  most  in- 
timately related  to  its  scattering  characteristics  at  low  frequencies. 
In  fact,  at  the  highest  frequency  of  interest  the  maximum  linear  di- 
mension of  the  target  is  only  one  and  one-half  to  two  wavelengths. 

That  is,  maximum  gross  physical  property  information  is  contained 
in  this  frequency  range.  From  the  viewpoint  of  ease  of  implementation 
obviously  higher  frequencies  are  of  interest.  One  must  recognize 
however  that  the  use  of  such  higher  frequencies  necessarily  entails 
some  loss  of  fundamental  identification  potential.  The  scattering 
characteristics  of  the  target  at  higher  frequencies  are  related  to 
isolated  features  of  the  target  rather  than  the  target  as  a whole. 

This  is  the  reason  why  approaches  to  the  target  identification  problem 
using  modulated  short  pulse  radars  have  generally  failed.  One  indeed 
isolates  the  flare  spots  by  this  approach  but  how  these  fit  together 
to  form  the  total  target  is  most  difficult  if  not  impossible  in  most 
cases  to  answer.  A bistatic  angle  for  the  radar  is  clearly  a designer 
controlled  parameter.  In  fact,  if  several  receivers  are  employed 
then  tracking  of  the  flare  spot  variations  with  target  orientation 
has  been  postulated  as  a modulated  short  pulse  radar  identification 
process.  The  intent  here  is  not  to  evaluate  this  or  any  other  iden- 
tification technique,  and  over  the  years  various  exotic  schemes  for 
coded  variations  of  the  controlable  parameters  have  been  proposed. 

We  note  however  that  as  the  class  of  targets  to  be  identified  grows 
larger  the  training  set  for  identification  by  this  approach  could 
quickly  become  unwieldy. 


It  is  customary  in  almost  any  diagnostic  or  identification  prob- 
lem to  utilize  all  of  the  available  observables.  The  approach  adopted 
here  however  minimizes  the  parameters  used  to  characterize  the  target. 
Each  target,  even  for  such  complicated  structures  as  aircraft,  is 
uniquely  characterized  by  a small  finite  stt  of  numbers.  These  num- 
bers are  the  difference  coefficients  of  a difference  equation  asso- 
ciated with  the  dominant  complex  natural  resonances  of  the  target 
or  (with  latest  research)  possibly  less  dominant  substructure  reso- 
nances. In  effect,  the  singularties  in  the  target  transfer  function 
are  modelled  as  simple  poles  and  are  excitation  invariant.  That 
is,  the  poles  are  independent  of  the  orientation  of  the  target  with 
respect  to  the  radar  and  of  the  polarization  of  the  interrogating 
radar  signal-3’  . It  follows  that  the  difference  equation  and  the 
difference  coefficients  have  the  same  properties. 

Processing  for  detection  or  identification  consists  of  comparing 
a measured  transient  response  from  an  unknown  target  with  a calculated 
transient  response  produced  using  the  coefficients  of  a known  target. 
Several  things  are  involved.  First,  the  measured  transient  response 
is  produced  synthetically  from  multiple  frequency  scattering  data 
which  with  substructures  can  perhaps  be  "chirp"  -type  signals.  Second, 
the  calculated  transient  response  is  obtained  using  a linear  difference 
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equation  and  the  coefficients  of  a given  target.  Thus  the  approach 
assumes  that  the  poles  or  coefficients  for  a class  of  target  shapes 
of  interest  have  previously  been  obtained.  The  fact  that  our  identi- 
fication process  does  not  extract  resonances  from  the  response  wave- 
form of  an  unknown  target  turns  out  to  be  a distinct  advantage.  The 
extraction  of  complex  natural  resonances  from  calculated  or  measured 
scattering  or  radiation  data  is  neither  a simple  nor  exact  task. 

The  resonances  can  be  found,  but  procedures  which  empirically  ranove 
anomalous  resonances  must  be  employed.  Data  from  a single  object 
orientation  and  a fixed  polarization  are  not  sufficient.  The  result 
is  that  resonance  extraction  is  not  a real  time  process.  Fortunately 


with  our  Identification  process  preparation  of  the  target  library 
is  a separate  task  and  not  directly  part  of  the  scheme.  Therefore 
our  method  does  yield  a real  time  capability.  Note  also  that  the 
basic  question  is  if  the  target  response  approximately  satisfies  a 
linear,  homogeneous  difference  equation. 

THE  BASICS 

The  impulse  response  concept  has  been  detailed  in  a tutorial 
2 

paper  . For  our  purpose  here,  it  suffices  to  note  that  the  transform 
pair 

G<M,p,ju)  « r Fr(9,<}>,S,t)e_jll)tdt 
0- 

A _ oo  A t 

Fr(9,<J),P,t)  = ^ / G(9,<M,jw)eJ+utdu,  m 


t 


relate  the  normalized  scattered  field,  G,  from  the  target  to  the  im 

A 

pulse  response  waveform,  Fj,  of  the  target.  The  triplet  9,  i>,  p 

denotes  the  spatial  and  polarization  dependence.  Our  main  interest 
in  the  past  was  in  the  ramp  response  waveform 


A 

FR(M,P,t) 


2tt 


(2) 


In  Equation  (2),  noting  that  G can  at  most  increase  as  jw  , it  is 
clear  that  the  ramp  response  waveform  emphasizes  the  low  frequency 
scattering  characteristics  of  the  target.  It  is  also  clear  from  (2) 
that  information  on  that  target's  gross  phyical  properties  is  con- 
tained in  the  ramp  response  waveform.  The  diagnostic  features  in- 
herent in  the  impulse  response  waveform  of  a target  should  not  be 
overlooked.  At  a fixed  orientation  and  polarizat ion,  one  real  time- 
dependent  waveform  sums  up  the  scattering  properties  of  the  object 
at  any  frequency.  Via  convolution,  the  response  to  any  arbitrary 
interrogating  signal  can  be  quickly  obtained.  All  of  the  features 
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illuminated  by  spectral  limited  asymptotic  and  moment  methods  are 
contained  in  this  waveform  and  can  be  separately  identified.  The 
proper  testing  procedure  for  any  postulated  identification  scheme 
is  via  the  canonical  (impulse,  step  and  ramp)  response  waveforms  of 
the  object.  Radar  imagery  features,  i.e.,  actually  producing  iso- 
metric three-dimensional  pictures  of  the  target  from  radar  data,  are 
most  basically  considered  in  the  time  domain.  We  will  elaborate  on 


For  all  but  the  most  perverse  target*  (those  with  jump  discontinuities 
in  the  ramp  response)  an  upper  summation  limit  of  eight  or  ten  provides 
a good  estimate  of  the  ramp  response  waveform.  Thus  discrete  scat- 
tering data  at  five  odd  harmonics,  (3),  or  ten  even  and  odd  harmonics 
(4),  suffice  to  determine  the  ramp  response  waveform.  Strictly  speak- 
ing, the  fundamental  should  come  from  Shannon's  sampling  theorem  but 
we  generally  will  not  know  the  effective  duration  of  the  response. 

The  response,  of  course,  decays  exponentially  but  we  take  the  effec- 
tive duration  of  the  waveform  to  be  the  time  to  decay  to  roughly  two 
percent  of  its  peak  magnitude.  A general  rule  of  thumb  is  that  at 
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the  fundamental,  the  maximum  linear  dimension  of  the  target  should 
be  approximately  one-fifth  of  a wavelength.  The  imagery  possibilities 
follow  directly  from  the  ramp  response  waveform.  It  has  been  shown^ 
that  the  magnitude  of  the  ramp  response  waveform  is  directly  propor- 
tional to  the  cross  sectional  area  of  the  target  along  the  line  of 
sight,  rigorously  to  the  shadow  boundary  and  approximately  beyond. 


Thus  image  production  consists  of  relating  one  or  more  cross  sectional 
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areas  or  profile  functions  to  an  object  geometry  . This  can  be  ex- 
pressed in  very  elegant  mathematical  terms  but  the  end  result  is  the 
same.  Note  particularly  that  there  is  a limitation  in  that  at  certain 
of  the  interrogating  frequencies  the  object  is  electrically  small. 
Higher  frequencies  can  obviously  sharpen  the  image  but  without  at 
least  synthetic  low  frequency  information  there  is  no  image.  We  stress 
this  limitation  because  difference  equation  identification,  which 
started  with  essentially  the  same  scattering  data  needed  for  the  ramp 
response,  does  appear  capable  of  using  somewhat  higher  frequencies 
alone,  i.e.,  without  the  use  of  very  low  frequency  information. 

The  term  "image"  as  used  here  refers  specifically  to  a three-dimen- 
sional estimate  of  the  target  geometry  obtained  from  an  interpretation 
of  cross  sectional  area  (transverse  to  the  radar  line  of  sight)  in- 
formation. The  cross  sectional  area  information  comes  directly  from 
the  forced  portion  of  the  target's  ramp  response  waveform.  It  follows 
that  aliasing  errors  caused  by  a lack  of  low  frequency  information 
negate  the  image  production  unless  special  procedures  are  postulated. 

In  any  event  the  fundamental  appears  vital. 


The  impulse,  step  and  ramp  response  waveforms  of  a conducting 
spherical  scatterer  for  backscatter  are  shown  in  Figure  1.  The  wave- 
forms shown  were  produced  synthetically  from  Mie  series  calculations 
using  respectively  6*,  40  and  10  harmonics  with  a fundamental  sphere 
circumference  of  0.?5  wavelengths.  We  repeat  these  waveforms  in  an 
attempt  to  make  this  report  as  self  contained  as  possible.  Earlier 
specular  and  creeping  wave  contributions  are  easily  discernable  in 
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the  Impulse  response  waveform  as  has  been  discussed  elsewhere'.  Note 
carefully  that  deformations  of  the  spherical  geometry  beyond  the  shadow 
boundary  can  only  alter  the  impulse  response  waveform  for  times  greater 
than  one  delay  time.  This  type  of  diagnostic  potential  has  been  little 
exploited  in  time  domain  studies.  The  ramp  response  waveform,  shown 
in  Figure  1 for  the  sphere,  has  been  studied  extensively  for  many 
targets  for  three  principal  reasons: 

1)  It  can  be  produced  synthetically  from  measurenents  at  a 
very  few  frequencies  therefore  making  ramp  response  radars 
feasible,  if  not  state  of  the  art. 

2)  The  forced  portion  of  the  response  is  directly  proportional 
to  the  target's  cross  sectional  area  and  therefore  radar 
imaging  via  ramp  response  radars  is  a demonstrated  capability* 


The  transient  portion  of  the  ramp  response  can  be  well  approxi- 
mated  by  a finite  exponential  sum",  meaning  that  excitation 
invariant  dominant  complex  natural  resonances  can  be  used 
to  characterize  the  target®’*®’^. 


It  is  interesting  to  observe  that  an  exponential  sum  approximation 
of  the  impulse  response  waveform  of  the  sphere  for  transit  times  greater 
than  one  would  appear  to  require  many  terms  to  approach  a reasonable 
estimate. 

LOCATION  OF  THE  TARGET  POLES 

First  efforts  at  effecting  target  classification  via  excitation 
invariant  parameters  viewed  the  dominant  complex  natural  resonance 
(poles)  as  the  fundamental  elements.  In  fact,  recent  studies  on  Pronv's 
method  and  its  variations  continue  to  view  the  poles  as  the  desired 
end  result.  It  now  appears  that  a more  meaningful  characterization 
should  perhaps  be  couched  in  terms  of  the  difference  coefficients 
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of  a linear  difference  equation.  The  reason  for  this  and  the  relation- 
ship between  the  poles  and  difference  coefficients  will  be  discussed 
later.  In  the  following  we  first  describe  the  methods  and  procedures 
for  locating  the  dominant  poles  of  a target.  Subsequently,  the  dif- 
ference coefficient  approach  is  discussed. 


If  one  makes  the  approximation  that 

- N a 6 t 

VM.p.t)  ■ l A (6,*,p)e  n , 


or  equivalently 

Q-i 

l Cr(0,*,pHjo))r 

I 0 (ju>)q 
q=0  q 

then  in  Equation  (6), 


s„  * -%  s 3 0 • ,s 

are  the  complex  natural  resonances  of  the  target  and  are  the  roots 
of 

? 0_(jo>)q  * 0 . (S 

q=0  q 

What  we  have  in  fact  done  is  to  assume  that  at  low  frequencies  the 
electromagnetic  scattering  problem,  truly  a problem  Involving  a dis- 
tributed parameter  system,  can  be  modelled  as  a lumped  circuit.  We 

note  that  the  concept  expressed  in  (6)  and  (7)  was  first  suggested 

2 

in  1Qfi5  , several  years  In  advance  of  the  so-called  singularity  ex- 
pansion method*.  We  stress  again  the  approximation  involvea. 


I 


'r 


For  spheres  and  infinite  circular  cylinders,  the  poles  can  be 
obtained  via  solutions  of  transcendental  equations,  given  by  Stratton7 
for  spherical  targets.  For  an  arbitrary  target  geometry  other  tech- 
niques must  be  utilized.  The  first  attempt  to  locate  the  dominant 
poles  of  an  arbitrary  structure  was  by  rational  function  fits  (7) 
of  calculated  or  measured  multiple  real -frequency  scattering  data. 

That  is,  a system  of  simultaenous  linear  equations 


G(9,<j>,p,ju).) 


1 + 


q=l 


,^i>( 


Q-i 

+ l C (0,<M)(jV.)r 

r*l  r 1 


(10) 


was  solved  for  the  D„  and  C coefficients  and  then  the  roots  of  the 

q r 

polynomials  were  extracted.  Equation  (1.0)  was  applied  to  various 

q 

target  orientations  to  test  the  excitation  invariance  of  the  poles' . 
It  was  found  that  while  the  poles  did  not  display  the  precise  invar- 
iance predicted,  the  results  were  sufficiently  successful  to  warrant 
further  study.  The  rational  approximant  approach  was  temporarily 
abandoned  with  the  advent  of  an  iterative  search  procedure  applied 
to  an  integral  equation  formulation  of  the  scattering  or  radiation 
problem^’®. 


Details  of  a reaction  integral  equation  formulation  for  the  radi- 

q 

ation  or  scattering  by  thin-wire  structures  have  been  given  . For 
arbitrary  wire  structures  it  is  somewhat  simpler  to  obtain  the  poles 
using  a radiation  approach.  It  should  be  understood  that  one  must 
assure  that  all  of  the  dominant  poles  of  the  structure  have  been  ex- 
cited. This  is  done  by  using  several  asymmetrically  located  gener- 
ators on  the  structure.  The  integral  equation  is  reduced  to  the  fi- 
nite matrix  equation 

[Z(s)] [I (s) J = [V ( s )]  . (11) 

Z ( s)  is  an  n x n open  circuit  impedance  matrix  whose  size  is  dictated 
by  the  number  of  wire  segments  used,  [I(s)]  and  [V(s)]  are  current 
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and  voltage  column  vectors  and  (s)  Is  complex  frequency.  Formally, 
the  poles  are  given  by  the  solutions  of 

l[Z(s)]|-0.  (U) 

Numerically  it  is  somewhat  simpler  to  find  the  peak  magnitude  of  the 
admittance  Y(s)  at  the  generator  terminals.  If  s^  is  a complex  fre- 
quency starting  point  then 

Z(s,) 

s = si ' mn — (13) 

S I s = s1 

is  applied  iteratively  (s  - s^  a small  step)  until  further  changes 
do  not  occur.  The  numerical  search  procedure  is  attractive  in  the 
sense  that  it  will,  when  applied  correctly,  always  lead  to  a set  of 
poles.  However,  it  has  some  serious  disadvantages.  Modelling  of 
the  target  structure  by  a wire  grid  configuration  can  become  quite 

Q 

elaborate  for  complex  structures  such  as  aircraft  . Calculations 
at  complex  frequencies  using  Richmond's  program*®  are  normally  quite 
efficient  but  the  iterative  search  procedure  is  horrendously  inef- 
ficient. If  more  than  just  the  dominant  poles  are  needed,  the  prob- 
lems are  magnified  because  more  wire  segments  are  needed  to  model 
the  target  at  higher  real  frequencies.  At  this  time  there  is  also 
an  upper  frequency  limit  based  on  the  electrical  size  of  the  struc- 
ture. 

A third  method  of  locating  the  poles  utilizes  what  is  known  as 
11  4 

Prony's  method  ’ . Prony's  method  linearizes  the  process  of  obtain- 
ing an  approximate  fit  of  a waveform  by  a finite  sum  of  exponentials 
{Equation  (5)),  and  is  adequately  explained  in  the  literature.  A 
less  restrictive  method,  termed  an  eigenanalysis  solution*9,  can  also 
be  applied.  Both  can  be  explained  as  an  N term  difference  equation 
fit  to  sampled  (At)  data  yielding 


The  Prony  and  eigenanalysis  methods  are  two  ways  for  finding  the  com- 
ponents of  [x].  With  the  Prony  method,  the  xQ  coefficient  is  set 
to  unity  and  the  squared  error  minimum  obtained  with  this  constraint. 

The  eigenanalysis  approach  chooses  the  components  of  [x]  as  those 
of  the  eigenvector  of  the  minimum  eigenvalue  of  the  symmetric  matrix 
| A | . The  two  methods  lead  to  different  squared  errors  , indicating 
that  the  constraint  imposed  in  the  Prony  approach  is  real,  i.e.,  some 

OC 

generality  is  lost.  It  has  recently  been  shown  that  there  are 
really  N+?  least  squared  error  solutions  for  the  coefficients  if  inter- 
polation methods  are  included. 

I 

S 

1 

The  Prony  and  eigenanalysis  methods  have  the  advantage  that  they 
can  be  applied  to  either  measured  or  calculated  transient  response 
waveforms.  Our  normal  procedure  has  been  to  obtain  a ramp  response 
waveform  via  Fourier  systhesis  using  (3)  or  (4)  and  either  calculated 
or  measured  data.  8oth  methods,  however,  do  have  some  problems. 

It  is  necessary  to  apply  them  over  different  later  (after  the  exci- 
tation Is  gone)  time  portions  of  the  waveform  until  a best  fit  is 
obtained.  It  does  not  necessarily  follow  that  the  poles  obtained 
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via  these  methods  will  be  the  same  as  those  obtained  using  the  in- 
tegral equation  approach.  An  example  of  this  will  be  given  in  another 
section  of  the  text.  Both  methods  yield  almost  amazingly  good  fits 
to  the  Intermediate  and  late  time  portions  of  transient  signals  but 
no  physics  is  associated  with  the  algorithm.  Thus,  occasionally, 
right  half-plane  poles  crop  up  as  do  poles  which  give  an  excellent 
fit  to  the  waveform  but  are  anomalous.  When  the  target  orientation 
and  radar  polarization  are  such  that  a pole  is  only  very  weakly  ex- 
cited, the  methods  seem  to  extract  the  imaginary  part  of  the  pole 

correctly  but  large  errors  can  occur  in  the  real  part.  This  also 

19 

occurs  when  the  data  are  noisy  . It  is  our  experience  that  best 
results  are  obtained  using  the  methods  when  all  of  the  desired  poles 
are  extracted  simultaneously.  Procedures  which  suggest  systematically 
working  toward  earlier  times,  i.e.,  increasing  the  time  span  covered 
and  correspondingly  increasing  the  number  of  poles  extracted  may  work 
well  with  ideal  data  but  are  not  effective  for  real  world  data  where 
the  errors  are  cumulative. 

Recently  the  rational  approximant  approach  to  target  pole  lo- 
cations has  been  applied  to  admittance  calculations  at  real  and  com- 
plex frequencies.  That  is,  the  region  of  the  complex  frequency  plane 
of  interest  is  spanned  by  a regular  grid  of  frequencies.  The  data 
were  then  utilized  as  described  by  (10)  and  the  subsequent  discussion. 

As  could  be  anticipated,  poles  extracted  via  this  approach  were  in 

very  good  agreement  with  those  previously  obtained  via  interactive 

1? 

search  or  the  Prony  and  eigenanalysis  methods  . The  improvement 
over  results  previously  obtained  using  only  real -frequency  data  was 
dramatic.  This  approach  to  target  pole  extraction  has  two  principal 
disadvantages.  First,  it  is  not  applicable  to  measured  data  and  second, 
for  a target  whose  poles  are  completely  unknown,  the  required  density 
of  complex  frequency  calculations  may  preclude  use  of  the  method  on 
complicated  target  structures.  In  any  event,  the  calculations  will 
have  to  be  computationally  efficient. 
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EXAMPLES  OF  DOMINANT  TARGET  POLES 


In  this  section  we  will  present  dominant  target  poles  for  a 
number  of  structures,  both  simple  and  complex.  Examples  will  be  given 
of  poles  found  by  certain  of  the  methods  discussed  in  the  previous 
section.  We  concern  ourselves  here  with  relationships  which  can  be 
found  between  the  target  size  and  shape  and  the  pole  locations.  The 
reader  should  note  that  in  every  case  only  the  second  quadrant  poles 
are  shown.  There  are  of  course  complex  conjugate  poles  in  the  third 
quadrant.  Also,  the  poles  are  plotted  in  the  complex  wavenumber  length 
plane  and  sometimes  with  an  additional  pi  normal ization. 

Figure  2 shows  the  first  four  dominant  pole  loci  of  a thin  wire 
(length  to  diameter  ratio  0.002)  as  it  is  bent  about  its  midpoint. 

The  resonance  at  lowest  frequency  in  Figure  2 is  identified  as  the 
first  even  resonance  because  at  this  frequency  the  current  distri- 
bution on  the  wire  is  even  with  respect  to  the  center  of  the  wire 
with  no  current  nulls  except  at  the  endpoints  of  the  wire.  For  the 
first  four  modes  the  current  distributions  resemble  respectively  one- 
half  cosine,  one  cycle  sine,  one  and  one-half  cycle  cosine  and  two 
cycle  sine.  With  this  in  mind  one  can  anticipate  that  as  the  wire 
is  bent  nearly  double  the  even  modes  will  migrate  to  the  imaginary 
axis  to  become  nonradiating  or  energy  storage  modes  because  they  can- 
not exist  on  a wire  one-half  the  length  of  the  original  wire.  This 
is  seen  in  the  loci  of  Figure  2.  Conversely,  the  odd  modes  will  mi- 
grate to  take  up  new  positions  as  respectively  first  even  and  first 
odd  modes  on  the  new  wire  length.  These  loci  terminations  are  shown 
as  crosses  in  Figure  2.  There  is  one  additional  lesson  to  be  learned 
from  Figure  2.  Proximity  of  the  pole  to  the  imaginary  axis  does  not 
necessarily  indicate  a strong  resonance,  one  must  also  know  the  resi- 
due of  the  pole. 


The  loci  shown  in  Figure  2 were  found  using  the  iterative  search 
procedure  with  one  or  more  generators  located  at  asymmetrical  points 
on  the  wire.  As  we  progress  to  more  complex  structures,  several  gen- 
erator locations  will  have  to  be  used  in  an  attempt  to  locate  al 1 
of  the  dominant  poles.  With  the  simplified  search  procedure  described 
in  the  previous  section  there  is  no  way  to  guarantee  that  all  the 
poles  have  been  found.  One  could  of  course  use  Cauchy's  residue  the- 
orem by  systematically  numerically  integrating  around  contours  in 
the  complex  plane  but  this  would  be  expensive  and  time  consuming  and 
the  search  procedure  is  already  inefficient.  As  will  be  seen  later, 
one  can  usually  detect  an  absent  pole  hy  applying  the  discrimination 
test  for  several  target  orientations. 

Figures  3 and  a show  the  poles  of  two  perpendicularly  crossed 
wires.  In  Figure  3 the  wires  are  crossed  at  their  centers  and  in 
Figure  A they  are  crossed  at  one-third  the  length  of  one  wire,  the 
pole  nearest  the  imaginary  axis  in  Fiqures  3 and  4 corresponds  to 
a bent  or  L-tyoe  current  mode  on  the  wires.  The  poles  in  Figures 
2 ( a = 90°),  3 and  4 show  a reassuring  self-consistency  in  that  the 
poles  corresponding  to  the  first  and  second  even  modes  are  undisturbed, 
as  should  be  the  case.  The  poles  in  Figures  3 and  4 were  also  found 
using  the  interative  search  procedure.  The  poles  of  two  simple  stick 
models  of  straight  and  swept  wing  aircraft  are  shown  respectively 
in  Figures  5 and  6.  Again,  the  interative  search  procedure  was  used 
to  find  the  poles  indicated  by  circles.  The  poles  shown  as  crosses 
will  be  discussed  shortly.  Note  here  however  that  only  one  pole  is 
significantly  different  for  the  two  aircraft.  We  shall  return  to 
this  ooint  later. 

The  dominant  poles  of  a conducting  circular  th’n  disk  of  radius 
a are  s^own  in  Figure  7.  These  poles  were  found  by  a method  suoqested 
hy  Berni  \ which  is  essentially  equivalent  to  Prony's  method.  The 
ooles  were  extracted  from  transient  waveforms  oroduced  synthetically 
from  calculated  frequency  domain  scattering  data'4.  The  results  in 
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Figure  7 show  one  basic  weakness  of  the  Prony  method.  The  poles  shown 
as  crossed  are  the  actual  poles  found.  However,  for  the  regular  disk 
geometry  one  would  anticipate  a string  of  poles  similar  to  those  for 
a wire  or  sphere.  The  imaginary  part  of  the  poles  show  this  type 
of  behavior  but  certain  of  the  real  parts  do  not.  The  proper  loca- 
tions for  these  poles  are  Indicated  by  circles.  That  these  positions 
were  correct  was  confirmed  by  comparing  transients  calculated  using  • 
the  poles  to  the  transients  obtained  synthetically  from  the  calculated 
frequency  domain  data.  It  has  been  our  experience  that  when  the  Prony 
method  falls  It  first  yields  incorrect  real  parts  for  the  extracted 
poles.  The  Prony  results  in  Figures  5 and  6,  which  differ  from  the 
interative  search  poles  except  for  one  pole,  were  not  the  result  of 
failure  of  the  Prony  method.  In  this  particular  case  the  frequency 
domain  scattering  data  used  to  synthesize  the  transient  response  only 
spanned  uL/ir  C = 0.1  ^0.1 11.0.  Clearly  the  frequencies  were  much 
too  low  to  excite  the  other  poles  found  via  interative  search.  This 
does  illustrate  another  type  of  difficulty  with  the  Prony  method. 

The  number  of  poles  to  be  extracted  must  be  specified,  and  it  is  not 
always  clear  from  the  residues  which  poles  are  anomalous. 

20 

The  dominant  pole  locations  for  eight  aircraft  have  been  given  . 

These  poles  were  obtained  via  Prony' s method  applied  to  transient 

scattered  waveforms  from  the  aircraft.  It  is  emphasized  that  the 

transient  waveforms  were  produced  from  good  electromagnetic  models 
g 

of  the  aircraft  . That  is,  calculated  backscatter  data  agreed  well 
with  experimental  model  measurements.  To  date,  most  of  the  identi- 
fication studies  using  predictor-correlator  processing  has  concen- 
trated on  the  F-4,  F-104  and  MIG-19  aircraft.  The  poles  for  these 
targets  are  given  in  Table  II. 


TABLE  II 


Aircraft 

o(xlO”6) 

jw  x 10"6 

F-104 

-19.7 

± 9 4.00 

- 7.01 

± 46.8 

- 3.5 

±131.0 

F-4 

-16.4 

± 66.0 

-10.9 

± 43.7 

- 9.9 

±130.0 

MIG-19 

- 8.7? 

± 79.3 

- 6.66 

± 55.6 

- 0.35 

± 66.7 

- 7.40 

±130.0 

Dominant  Natural  Resonances  of  F-104,  F-4 
and  MIG-1.9  Aircraft. 


It  should  be  noted  in  passing  that  once  the  complex  natural  reso- 
nances of  a class  of  targets  are  known,  there  are  a number  of  ways 
one  might  postulate  exploiting  them  for  identification.  Note  however 
that  with  present  resonance  extraction  methods  it  is  not  feasible 
to  suggest  real  time  extraction  of  the  resonances  from  a signal  ob- 
tained from  an  unknown  target.  Using  the  impulse  response,  the  re- 
sponse to  an  arbitrary  interrogating  signal,  E’(t)  is 


FSft)  = / F.ft-  rlE’fx  )dr 

o 


?A 


A method  of  rapidly  calculating  this  convolution  has  been  given^. 

To  visualize  the  possibilities,  approximate  the  impulse  response  as 
an  Infinite  summation,  then 


FS(t)  = l A 


n=l 


Je 
o 


Yn(t-r)  i 

n eV  x)dt 
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If  certain  of  the  resonances  (y  1 are  known  then  various  forms  of 

'n 

the  interrogating  signal  might  be  used  to  exploit  this  knowledge. 

This  is  perhaps  best  seen  in  the  frequency  domain  recalling  that  only 
a finite  number  of  poles  need  be  excited 


/rFS't1]  = l -A-  /[ fcVt)]  Mol 

n=l  Tn 

SUBSTRUCTURE  COMPLEX  NATURAL  RESONANCES 

It  has  long  been  recognized  that  a major  problem  associated  with 
the  successfully  demonstrated  dominant  complex  natural  resonance  ap- 
proach to  identification  developed  at  this  laboratory  is  the  spectral 
radar  data  required.  The  requirement  of  multiple  frequency  complex 
coherent  radar  data,  while  severe,  is  within  the  state  of  the  art. 
Moreover,  it  has  been  shown  that  certain  of  the  frequencies  and  phases 
can  be  deleted  . However,  the  fact  that  the  required  sampling  fre- 
ouencies  are  low,  e.g.,  roughly  2.0  to  20.0  MHz  for  fighter-type  air- 
craft, is  much  more  difficult  to  circumvent.  There  are  obvious  beam- 
width  and  other  problems  although  it  should  be  noted  that  if  a long 
ranqe  capability  were  sought  these  frequencies  for  an  over  the  horizon 
svstem  are  not  unrealistic.  In  this  case  there  are  other  problems 
such  as  coherency  for  the  bounce  path  which  would  require  study.  In 

anv  event  it  is  doubtful  that  an  airbourne  identification  radar  which 
exploits  dominant  complex  natural  resonances  of  targets  is  feasible 
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except  on  large  airframes  unless  a synthetic  aperture  approach  were 
employed.  For  this  reason  methods  have  been  sought  whereby  the  ex- 
ploitation of  target  complex  natural  resonances  for  identification 
is  retained  but  the  required  interrogation  frequencies  are  raised 
significantly. 

Any  finite  three-dimensional  object  has  an  infinite  number  of 
singularities  In  its  transfer  function,  l.e.,  that  function  which 
transforms  the  incident  field  into  the  scattered  field  at  a great 
distance  from  the  object.  Thus  in  principle  there  are  an  infinite 
number  of  complex  natural  resonances  (it  does  not  necessarily  follow 
that  all  of  these  singularities  can  be  modeled  as  poles  of  some  finite 
order).  These  resonances  are  basically  of  two  types  if  we  view  the 
second  quadrant  of  the  complex  frequency  plane  as  horizontal  1 avers 
bounded  by  zero  and  some  constant  values  of  oscillatory  frequency. 
There  are  1)  multiples  of  the  dominant  complex  natural  resonances 
and,  as  the  upper  bound  (nonzero  constant  jw)  increases,  2)  the  be- 
ginning of  new  strings  of  resonances  as  smaller  geometrical  features 
of  the  object  reach  complex  resonance  proportions.  If  the  layers 
(two  constant  values  of  jw)  are  thought  of  as  strips  of  constant  width 
then  it  follows  that  the  density  of  resonance  locations  within  a strip 
increases  as  ju>  increases.  A primary  objective  of  this  research  pro- 
gram was  to  determine  if  target  identification  could  be  achieved  using 
complex  natural  resonances  associated  with  smaller  geometric  features 
of  the  object.  This  choice  as  opposed  to  multiple  resonances  was 
made  because  it  is  known  that  multiples  of  a dominant  resonance  have 
an  ever-increasing  damping  with  respect  to  multiplicity,  i.e.,  the 
resonances  become  like  low-Q  structures. 

For  aircraft  the  two  most  logical  substructure  resonance  features 
are  combinations  of  the  vertical  and  horizontal  tail  stabilizers  and 
the  intake  cavities  of  the  jet  engines.  In  this  research  program 
both  features  have  been  studied.  The  approach  is  experimental  because 
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detailed  computations  which  must  necessarily  Include  the  Influence 
of  the  supporting  airframe  would  be  a study  In  itself.  Even  experi- 
mentally the  study  is  complicated  by  the  fact  that  for  these  features 
estimates  of  the  appropriate  spectral  region  where  the  feature  will 
be  most  influential  are  most  difficult.  This  Is  particularly  true 
for  the  stabilizers  since  both  (horizontal  and  vertical)  must  be  con- 
sidered as  supporting  certain  basic  modes.  The  philosophy  behind 
the  measurement  program  and  the  results  of  that  measurement  program 
are  given  in  a separate  report  . For  six  1/7?  scale  models  of  fighter 
aircraft  and  one  pencil-type  model  containing  a vertical  staMlizer 
and  a shorted  cavity,  swept  frequency  amplitude  and  phase  measure- 
ments were  made  In  appropriate  spectral  ranges.  From  these  data  syn- 
thetic pulse  response  waveforms  were  generated.  From  these  results 
it  is  possible  to  conclude  that  for  certain  aircraft  both  the  sug- 
gested substructure  features  have  at  least  partially  dominating  in- 
fluence in  the  approximate  full  scale  range  of  28.0  to  TOO  MHz.  How- 
ever, a high-Q  type  resonance  effect  was  not  observed.  The  processing 
of  these  data  and  the  results  are  reported  here.  The  basic  question 
now  being  if  an  excitation  invariant  linear,  homogeneous  difference 
equation  can  be  found  for  at  least  one  target,  the  implication  being 
others  can  be  similarly  treated. 

A final  point  needs  to  be  made,  particularly  since  only  very 
limited  results  are  reported.  The  basic  interrogating  frequency  range 
for  identification  continues  to  be  that  span  where  the  maximum  linear 
dimension  of  the  target  is  at  most  a wavelength.  !t  is  this  span 
of  frequencies  which  elicit  the  fundamental  physical  properties  of 
the  ta^aet.  It  follows  therefore  that  if  the  interrogatina  frequencies 
are  raised  above  the  fundamental  span  as  is  suggested  here,  one  can 
anticipate  a possible  corresponding  degradation  in  target  identification 
results.  Identification  results  using  measured  noisy  data  in  the 
fundamental  range  were  excellent4.  Therefore  good  identification 
may  still  be  possible  using  substructure  features  but  some  possible 
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loss  of  potential  may  ensue.  As  noted  earlier,  it  was  found  experl- 
mentally  that  a response  from  the  tail  region  dominated  the  return 
from  the  aircraft  in  the  30.0  to  56.0  MHz  region  but  the  characteris- 
tics of  a hlgh-Q  resonance  were  not  observed.  In  many  cases  only 
a single  reflected  pulse  is  seen.  We  seek  here  to  use  the  envelope 
of  this  pulse  response  beyond  the  peak  of  the  pulse  to  characterize 
the  target  via  a difference  equation.  It  follows  that  "poles"  ob- 
tained from  the  resulting  difference  coefficients  will  not  have  a 
physical  interpretation  in  terms  of  the  stabilizers.  This,  however, 
does  not  preclude  using  the  coefficients  for  Identification. 

PROCESSING* 

With  the  type  of  analysis  involved,  the  basis  question  is  if 
a portion  of  a discretized  transient  record  approximately  satisfies 
a linear  homogeneous  difference  equation 

N 

n=0  ^N+m-nan  ~ em’ 

; ms0,l,2, . . .M;  M+l  > N+l  . (201 


where 

fi  - ffiAt)  1211 

are  the  data  samples  taken  at  the  sample  interval  At.  The  t are 

m 

t.He  errors,  assumed  small,  and  the  ap  are  the  unknown  coefficients. 

The  equations  in  (’0)  can  be  put  in  matrix  form  as 

fF)(A!  = f e) , 1221 


♦Certain  of  the  ideas  summarized  here  were  obtained  from  research 
on  another  program^. 


where  (F)  is  an  (M+l)x(N+l)  data  matrix,  (A)  an  N+l  component  column 
vector  and  (e)  an  M+1  component  error  vector.  The  total  squared  error 
is 

Ej  * U)T(e)  (23) 

where  (e)T  denotes  transpose.  The  total  squared  error  is 

Ej  * (A)T(F)T(F)(A),  (24) 


where  the  right  side  of  (24)  is  a standard  quadratic  form.  Define 
(R)  = (F)V),  (25) 

then 

Ej  - fA)T(R)(A),  (26) 

where  (R)  is  a real  symmetric  (N+l)x(N+l)  matrix.  There  are  N+2  dif- 
ferent minimum  total  squared  error  solutions  for  the  coefficients. 

The  most  qeneral  in  that  it  imposed  no  restrictions  on  any  one  of 

iq 

the  coefficients  is  that  of  eigenanalysis  . Since  fRl  is  real  and 
symmetric 

eJ)  = (AlT(P)T(A)(P)(Al,  r 27) 

where  (A)  is  a diagonal  matrix  of  the  eigenvalues  of  (A).  Then  with 
the  transformation 

(Y)  - (P) (A) , (281 

and 


With  the  restriction 


||A||-1,  (30) 

The  eigenvalues,  N+l  in  number,  are  minimum  squared  errors  and  the 
minimum  eigenvalue  is  the  minimum  squared  error.  The  remaining  eigen- 
values are  minimum  squared  errors  with  additional  fin  addition  to 
(30))  constraints  in  terms  of  inner  products  of  the  associated  eigen- 
vectors. The  eigenanalysis  approach  was  also  formulated  in  the  z- 

transform  domain.  Based  on  this  a general  optimization  procedure 

19 

was  tested  on  another  program  . Despite  the  fact  that  a z domain 
approach  seems  ideal  for  discrete  data,  no  particular  advantage  could 
be  found  with  this  method.  Accordingly,  further  study  of  z-domain 
methods  are  abandoned.  The  roots  of  the  characteristic  equation  ob- 
tained from  the  coefficients  are  excitation  invariant  and  it  follows 
that  the  coefficients  are  excitation  invariant.  This  suggests  that 
data  records  corresponding  to  several  aspects  and/or  polarizations 
for  a qiven  target  can  be  combined.  Assume  that  we  have 

0 = \?,...0,  f3\1 

tarqets  and  for  each  target  we  have 

1 * 1,2,.. .L*,  (32) 

records.  For  each  record  and  target  a real  symmetric  matrix 


(R<*>)  . I F<l>)T 


(33) 


♦Clearly  in  general  L=Lq  and  M=M^. 


1 in  wiiwjyprin  nu 


Is  obtained.  Now  define 


% > '%  > • 


and  (Rq)  can  be  used  for  (R)  In  (26).  Since  the  individual  records 
can  have  different  effective  S/N  ratios  a weighting  of  the  R^  is 
also  possible.  This  modified  elgenanalysis  approach  has  been  used 
for  the  data  obtained  on  this  contract. 


The  remaining  N+l  methods  for  finding  the  coefficients  consist 
of  normalizing  (26)  by  one  of  the  coefficients 


rf  fA)T  (R) (A)  , 


v 


where  v is  an  integer  (see  (’01)  and 


0 < v < N , 


E?/av 

6 4-^  = 0,  J7v 

to. 

is  used  to  obtain  the  normal  equations.  For  example,  for  Prony's 
method  (20)  becomes 

N a„  £m 

c A r f JI 


+ y f JL  = 

nti  N+m-n  a„  art  * 

n=i  o o 


m=0,l,...M 


where  the  normalization  by  aQ  is  usually  absorbed. 


In  general 


n cJ  _ 

f + r r J2  - em 
N+m-v  N+m-n  a..  ~ a 


m=0,t,?,...M. 


For  proper  v,  the  squared  error  obtained  with  Prony's  method  can  be 
reduced  by  almost  a factor  of  ten.  A general  v approach  Is  being 
studied  on  another  program  . The  general  v approach  was  suggested 
however  because  of  Identification  using  prediction-correlation.  As- 
sume Prony's  method  is  used  to  find  the  coefficients  and  one  has 


^N+m  + ^N+m-n  an  ~ 


(40) 


where  the  normalization  has  been  absorbed  and  zero  error  now  assumed. 
With  the  coefficients  known,  (40)  is  a predictor,  i.e.,  fN+m  can  be 
predicted  from  N previous  samples.  But  If  the  samples  are  noisy  then 
(40)  should  not  be  used  as  is  because  low  S/N  samples  can  be  multiplied 


by  large  numbers  (coefficients), 
ficient  then  ^40)  is  normalized  as 


Instead  if  ay  is  the  largest  coef- 


N+m  f N a 

~ N+m-y  + T fw  — = 0. 
ay  n4i  N+m-n  - 

ni'y  y 

The  predictor  equation  becomes 


(ail 


N+m 


N+m-y 


a ' n«l  f*"-»  ^ • 
y nty  ay 


(42) 


Now  note  that  in  general  (42)  is  an  interpolation,  i.e.,  fN+n  comes 
from  both  prior  and  future  samples.  The  equivalent  expression  if 
the  coefficients  are  found  by  eigenanalysis  is 


N-y 


N 

I '■ 

n=0 
n/y 


N-n  . 

y 


(43) 


In  both  (42)  and  (43),  all  of  the  coefficients  multiplying  noisy  data 
have  been  reduced  in  magnitude.  The  next  section  of  this  report  dis- 
cusses the  general  Identification  problem. 
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IDENTIFICATION 


The  target  characterization  parameters  are  the  excitation  in- 
variant coefficients  of  the  linear  homogeneous  difference  equation 
which  is  satisfied  by  scattered  data  from  the  target.  Equivalently 
one  could  say  that  the  characterization  parameters  are  the  excitation 
invariant  complex  natural  resonances  of  the  target  or  that  the  target 
is  characterized  by  the  excitation  invariant  difference  equation. 

The  tarqet  is  also  characterized  by  the  excitation  dependent  residues 
in  (6)  and  (71.  Use  of  the  residues  however  would  require  much  larger 
storage  libraries  and  a pattern  recognition-type  learning  procedure. 

We  have  consistently  advocated  a procedure  which  does  not  require 
parameter  extraction  from  an  unknown  (sourcel  waveform  resDonse. 


Let  f;*'  be  samples  of  the  response  from  an  unknown  target  and 

( c ^ 

fj  be  a calculated  waveform.  First,  the  calculated  waveform  is 
obtained  using  (43)  which  for  general  m is 


f(c)  = _ y f(t)  ji 
N+m-y  N+m-n  a * 

nj<y 

m=0,l,...M. 


(44) 


Thus  a calculated  waveform  is  obtained  using  the  coefficients  of  a 
known  waveform  (from  the  characterization  library)  and  samples  of 
an  unknown  waveform.  This  is  prediction.  One  can  now  calculate 


x » 1 - 


M 


(t) 

N+i 




y tf 

1-0 


f th2 
N+i' 


p(c) 
rN+i ' 


(f 


fc) 

N+i 


2 


(45) 


which  is  unity  minus  the  normalized  squared  error  between  the  test 
and  calculated  waveforms  for  the  fixed  value  of  the  samplinq  interval. 
The  characteristic  equation  associated  with  (?0)  is 


? VN'"  = 0.  (46) 

n*0 

the  roots  of  which  are 
y.At 

s.  = e 1 , 1-0,... N.  (47) 

Thus  the  complex  natural  resonances  are 

Yl’Hlnfs1)*  m) 

Equatinq  f46)  and  a factored  form  of  (46) 


N 

1 V 

n=0  n 


N-n 


YjAt 

= ao  i3l  'S+e  > ’ 


N 


faq) 


yields  how  the  coefficients  must  change  with  At.  Finally 


P"  * I 1 - x(At) , f?0) 

At 

is  the  identification  function.  If  the  test  waveform  is  from  the  target 
whose  coefficients  have  been  used  then  values. of  p"  near  unity  are 
obtained.  Low  values  of  p"  indicate  a mismatch  between  the  test  wave- 
form and  target  coefficients. 


There  are  of  course  a number  of  other  methods  whereby  the  pa- 
rameters discussed  above  might  be  used  in  identification.  However, 
the  prediction-correlation  shown  has  performed  satisfactori Iv.  For 
our  present  purpose  the  major  question  is  the  feasibility  of  obtaining 
reasonably  stable  parameters  usinq  data  associated  with  substructure 
features  but  not  resonances  of  the  target.  The  actual  processing 
results  obtained  for  pulse  response  of  aircraft  via  (14)  are  discussed 
in  a later  section  of  this  report. 
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EIGENANALYSIS  APPLIED  TO  MEASURED  DATA 

Another  report  on  this  contract-^  describes  methods  of  obtaining 
measured  spectral  and  temporal  scattering  data  in  some  detail.  Bas- 
ically, sampled  frequency  domain  data  (magnitude  and  phase)  were  meas- 
ured in  two  spectral  regions;  2. 2-4.0  GHz,  and  4. 0-7. 6 GHz  for  several 
aircraft  models  and  a pencil-type  cavity  structure.  Each  spectral 
band  was  spanned  by  201  evenly  spaced  points  (200  intervals).  Although 
data  from  both  frequency  bands  were  reported  on  Reference  24,  eiqen- 
analysis  was  only  performed  on  the  2. 2-4.0  GHz  data.  Time  did  not 
permit  more  than  a cursory  examination  of  the  higher  frequency  results. 

In  previous  work  on  identification  at  this  Laboratory,  transient 

time  waveforms  were  simply  obtained  from  ten  frequency  harmonically 

related  data.  Using  higher  and  more  realistic  frequencies  on  this 

contract  necessitated  a new  technique  for  obtaining  a suitable  time 

domain  waveform.  The  approach  employed  was  to  consider  the  2. 2-4.0 

GHz  data  to  be  modulated,  i.e.,  effective  demodulation  of  the  data 

from  -1.8  to  1.8  GHz  and  conversion  to  the  time  domain  via  Fourier 

Series  yielded  what  was  considered  to  be  the  "envelope"  of  a short 

pulse  return.  These  time  domain  pulse  returns  are  given  in  the  data 
74 

report'  and  discussed  there.  In  this  report  the  time  responses  are 
considered  as  the  data  bank  to  which  the  eiqenanalysis  method  was 
applied. 

One  additional  note  on  the  data  processing  is  made.  After  the 
frequency  domain  data  were  shifted  down  to  be  centered  at  zero  fre- 
quency a filter,  specifically  a cos  n window,  was  applied. 


The  cos  a window  can  be  expressed  as 


23 

The  value  of  a in  (51)  can  be  an  integer  , but  in  our  application  the 
value  of  a * 25/46  was  used. 


Method 

Equation  (?0)  given  earlier  views  the  general  process  in  terms 
of  the  difference  equation  for  a constant  sample  interval.  More  de- 


tail must  be  given  at  this  time  as  a slightly  different  procedure 


was  used  to  form  the  data  matrix  referred  to  as  (F)  in  Equation  (22). 
Rewriting  (20)  as  a function  of  time 

l0  "V"  ‘-n  l>an  - cm 

; m=0,l,...M;  (M+l)  > (N+l). 


It  is  seen  that  as  long  as  the  time  increment.  At,  multiplying  n in 
(52)  is  constant  in  all  equations,  the  difference  equation  method 
holds.  The  At  multiplying  m need  not  be  equal  to  the  At  multiplying 
n.  The  situation  can  be  seen  in  Figure  8.  Forcing  the  two  At's  to 
be  equal  is  the  situation  depicted  in  (8a);  no  overlapping  occurs. 

In  fRb),  however,  the  At  multiplying  m is  smaller,  the  At's  overlaD, 
and  more  use  can  be  made  of  the  appropriate  part  of  the  waveform. 
Calling  the  smallest  increment  between  samples  At,  Equation  (52)  can 
He  rewritten  as 

N 

f(tN+m6t-nAt)an  = em.  '53) 

Obviously, 


At  * L’St;  l an  integer  > 0. 

Then 
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-1  r8t 
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-i  h»»  / 
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*•  = At  = 3 8t 

OVERLAPPING  At's 


(b) 

Figure  8.  Nonoverlapping  (a)  and  overlapping  (b)  subintervals. 
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! ~ Any  given  value  of  m still  gives  one  difference  equation  formulation 

with  a constant  At;  the  At's  in  all  equations  are  equal.  We  have 
the  same  situation  in  Equation  (54)  as  in  (5?)  but,  again,  the  critical 
portion  of  the  waveform  is  possibly  employed  more  efficiently. 

~ Rewriting  (54)  as  in  Equation  (?0)  we  have 

| N 

n=0  fN+m"nL  an  = em* 

m = 0,1,?, ...,M;  fM+l)>(N+l). 

_ L an  integer  >0  (55) 

where  the  time  increment  is  now  6t. 

The  formulation  of  the  data  matrix  (F)  as  given  by  Equation  f?0) 
is  (allowing  At  to  be  equal  to  L) 

N 

n=0  ^N+mL-nL  an  em»  <56) 


and 


The  formulation  for  Equation  (551  Is 
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fN-L 

fN-2L 

f N-NL 

fN+l 
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(FR) 

• 

fN+L 
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• 

fN 

fN-L 

•”  fN+L-NL 

(L+11+h  row 

• 

fN+M 

fN+M-L 

fN+M-2l 

f N+M-NL 

These  two  formulations  are  not  simply  expressed  but  the  important 
feature  is  that  the  (L+l)  row  in  the  second  method  is  the  same  as 
the  2nd  row  in  the  first  method.  Again  it  is  seen  that  more  informa- 
tion can  possibly  be  obtained  from  the  important  portion  of  the  wave- 
form. 

The  new  procedure  exploited  in  this  study  was  the  use  of  infor- 
mation from  two  data  records  (i.e.,  two  aspect  angles!  in  forming 
the  data  matrix  (FI.  In  practice,  (F)  was  actually  formed  from  two 

submatrices.  Consider  one  target  and  a time  domain  curve  FVtl  from 

1 ? 

each  of  two  aspect  angles.  Call  them  F and  F . Several  constraints 
exist  in  forming  the  new  data  matrix  namely;  At  must  be  the  same  for 
both  curves,  the  number  of  terms  in  the  desired  difference  equation 
must  be  the  same  for  both  curves  and  the  number  of  rows  of  (F)  must 
be  greater  than  or  equal  to  the  number  of  columns.  This  last  constraint 
reflects  the  least  square  procedure  in  Equations  (23)  through  (30) 
with  reqard  to  the  (M+l)  _>  (N+l)  condition  of  Equations  (20)  and  (55). 
An  N term  difference  equation  (F)  could  be  composed  of  Ml  rows  of 

. 1 o 

information  from  F and  M2  rows  from  F . As  stated  above 


w+n+(M?+i.)  > n+i  rs9) 

but  note  that  Ml  is  not  necessarily  equal  to  M2.  We  now  represent 
the  data  matrix  for  the  following  difference  equation  system 
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The  formulation  given  by  (60)  is  for  two  data  records.  In  fact, 
the  method  could  be  applied  to  an  arbitrary  number  of  data  records. 

The  examples  given  in  this  report  were  restricted  to  two  aspect  angles. 


In  applying  Equation  (60)  to  our  measured  radar  data,  several 
parameters  had  to  be  varied.  They  are: 
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1. 


N,  the  number  of  terms  in  the  difference  equation. 


2.  L,  the  At  in  the  difference  equation. 

3.  Ml,  the  number  of  rows  of  information  from  data  record  F^. 

4.  M2,  the  number  of  rows  of  information  from  data  record  F . 

In  general,  another  variable  exists.  The  portion  of  the  wave- 
form over  which  the  data  is  sampled  can  be  changed.  In  this  study, 
this  variable  was  in  effect  removed.  Examination  of  (58)  reveals 
that  the  sample  corresponding  to  the  earliest  point  in  time  is  in 
the  upper  right  hand  corner,  FN_^.  the  synthesized  pulse 
responses  to  which  this  method  was  applied,  F^_^  was  always  considered 
to  be  the  maximum  point  on  the  waveform.  (See  discussion  on  pulse 
responses  in  section  on  substructure  resonances.)  The  sampled  "win- 
dow" of  the  waveform  was  always  the  maximum  point  and  points  to  the 
right  of  it  corresponding  to  longer  times.  Some  variation  in  the 
sampled  "window"  occurred  as  N and  L (or  At)  were  changed.  As  L was 
increased  for  a given  N,  samples  corresponding  to  longer  times  were 
included.  Increasing  N,  the  number  of  terms  in  the  difference  eq- 
uation, also  forced  inclusion  of  more  data  points.  Increasing  M, 
the  number  of  rows  in  the  data  matrix  corresponding  to  a given  record 
also  changes  the  size  of  the  "window".  The  starting  point,  however, 
was  always  the  maximum  point  on  the  waveform.  It  was  always  included 
in  this  study. 

As  these  variables  were  adjusted  the  characteristic  being  observed 
was  the  stability  of  the  difference  equation  coefficients,  the  a^, 
in  Equation  (601.  By  stability  here  we  mean  invariance  to  changes. 

In  this  application  for  a given  number  of  terms  in  the  difference 
eauation  N,  and  a fixed  At,  Ml  and  M2  were  varied  and  the  coefficients 
examined  for  stability.  Obviously,  a change  in  At  or  N will  yield 


totally  different  sets  of  coefficients.  The  actual  procedure  used 
then  was  to  pick  a value  of  N,  pick  a value  of  At,  and  allow  Ml  and  M2 
to  vary. 

Results  of  this  application  to  several  waveforms  are  given  in 
the  next  section. 

Results 

As  stated  earlier,  the  goal  of  this  portion  of  the  effort  was 
considered  to  be  stability  of  the  coefficients  for  a given  number 
of  terms  in  the  difference  equation  and  a given  At  as  more  information 
was  added  to  the  data  matrix  (F).  Again,  changing  information  in 
(F)  corresponds  to  adding  rows  of  elements  as  defined  by  a given  value 
of  m in  Equation  (55)  for  either  of  the  two  aspect  angles.  Equation 
(60)  defines  the  process  for  two  aspect  angles.  It  is  seen  that  there 
are  MM1  rows  of  data  from  aspect  angle  1 and  MM2  rows  of  data  from 
aspect  angle  2.  The  following  discussion  considers  NN  to  te  the  number 
of  terms  in  the  difference  equation  and  IDT  to  be  At.  The  "error" 
discussed  refers  to  the  total  squared  error  as  defined  by  (20),  (i.e., 
the  smallest  eigenvalue),  divided  by  the  trace  of  (R)  where  (R)  was 
defined  in  Equation  (25),  and  the  trace  of  this  real  symmetric  matrix 
is  the  sum  of  the  diagonal  elements.  This  weighting  permits  a fair 
comparison  between  errors  of  different  size  matrices. 

Application  of  this  procedure  proved  time  consuming  because  a 
systemic  procedure  was  not  developed  for  this  test.  Searchinq  for 
stable  coefficients  initially  involved  lengthy  computer  printouts 
and  visual  examination.  Figure  9 shows  a semi-automated  method  where 
the  coefficient  values  are  plotted  as  functions  of  their  coefficient 
numbers.  In  fact,  the  computer  program  has  the  capability  of  pro- 
ducing these  plots  on  a CRT  screen.  This  greatly  aids  examination 
of  the  output  data. 


Figure  9 shows  difference  equation  coefficients  plotted  when 
the  procedure  was  applied  to  the  time  waveforms  of  Figure  10.  (A 
pencil  shape  target  with  a vertical  tail  stabilizer  as  described  in 
Reference  24.)  As  shown,  seven  coefficients  were  sought  with  IDT=4. 

Note  that  as  MM1  and  MM2  are  varied  (the  size  of  the  data  matrix  changed), 
the  coefficients  vary  somewhat,  i.e.,  are  not  stable.  Actually,  vari- 
ations were  larger  than  this  plot  indicates.  Figure  R shows  only 
those  cases  where  the  first  coefficient  was  positive.  Figure  1.1  shows 
the  coefficients  plotted  with  six  terms  in  the  difference  equation 
and  IDT=4  for  the  same  pulse  response.  Again,  stability  is  not  in- 
dicated. The  case  of  seven  terms  in  the  difference  equation  and  IDT=3 
is  shown  in  Figure  12.  For  these  ten  values  of  MM1  and  MM2  stability 
is  clearly  achieved.  In  fact,  there  were  more  than  ten  cases  which 
would  have  shown  this  stability.  Figure  12  also  reveals,  however, 
that  the  first  coefficient  is  very  near  zero.  This  indicated  that 
perhaps  only  six  terms  are  actually  needed  in  the  difference  equation. 
Indeed  Figure  13  for  NN=6  and  IDT=3  indicates  stability  comparable 
to  Figure  12  with  no  near-zero  term.  The  resulting  difference  eq- 
uation indicated  by  Figure  13  which  best  fit  the  two  curves  of  Fig- 
ure 10  is 

.182  F V t- A)+ .0975F 1 ( t-2 At! - .0255F  * ( t-3At! 

+ .495F1(t-4At)-.722Fi(t-5At)+.354F1(t-6Atl  = 0 («2) 

where 

At  = 3*5  = 3x0.139  = .417  ns 
.139  ns  * spacing  between  time  samples 


1*1,2  refers  to  the  two  data  sets. 
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Figure  13.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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Figure  12.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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Figure  10.  Pulse  responses  of  pencil  target  2 inch  tail  stabilizer. 
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Figure  9.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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This  procedure  has  been  applied  to  several  other  waveforms  but  it 
is  interesting  here  to  examine  the  waveforms  of  Figure  14.  These 
waveforms  correspond  to  the  same  pent  shaped  target  on  Figure  10  with 
the  exception  that  the  vertical  tail  fin  is  slightly  longer.  Figure 
15  shows  a stable  set  of  coefficients  with  five  terms  and  IDT=4.  Fig- 
ure 16  shows  a slightly  better  fit  with  regard  to  the  total  squared 
error.  The  difference  equation  for  the  target  of  Figure  14  would 
be 

.718F1(t-At)-.289F1(t-2At)-.186F1(t-3At) 

-.451F1 ( t-4At)+.404Fi ( t-5At)  = 0 (63) 

These  two  examples  are  used  ( compare  Equations  (62)  and  (63) 
to  show  that  two  very  similar  targets,  differing  only  in  a single 
substructure  dimension,  yield  two  vastly  differing,  yet  stable,  dif- 
ference equations.  Past  experience  with  a pole  or  difference  coef- 
ficient discriminator  as  in  Equation  (45)^’^  indicates  that  discrimi- 
nation can  be  achieved  using  these  coefficients. 

DISCUS ' ! ON  AND  CONCLUSIONS 

The  early  sections  of  this  report  were  devoted  to  a summary  of 
a basic  method  of  target  identification  developed  over  the  course 
of  several  contracts.  These  developments  are  translated  here  to  a 
set  of  conclusions  addressed  specifically  to  modern  fighter  aircraft. 

1.  Scattering  data  in  the  spectral  range  of  2.0  to  24.0  MHz 
contain  the  most  basic  information  concerning  the  overall 
size  and  shape  of  the  aircraft. 
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Figure  14.  Pulse  responses  of  pencil  target  2.5  inch  tail  stabilizer. 
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Figure  15.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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Figure  16.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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Figure  17.  Difference  equation  coefficients  as  a function  of 
equations  from  two  aspects. 
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Transient  response  waveforms  can  be  produced  from  ten  dis- 
crete frequency  samples  of  the  complex  echo  siqnal  in  the 
2.0  to  24.0  MHz  range  which  are  relatively  simple  in  shaoe 
and  yet  characteristic  of  the  aircraft. 

The  transient  response  waveforms  produced  from  samples  in 
the  2.0  to  24.0  MHz  frequency  range  display  substantial 
ringing  indicative  of  the  presence  of  one  or  more  high-Q 
complex  natural  resonances. 

Excitation  invariant  complex  natural  resonances  can  be  ex- 
tracted from  transient  response  waveforms  produced  from 
complex  scattering  data  in  the  frequency  range  2.0  to  24.0 
MHz.  These  dominant  complex  natural  resonances  are  finite 
in  number  (usually  no  more  than  five  of  six)  and  are  phys- 
ically related  to  large  structures  such  as  the  wings  and 
fuselage  of  the  aircraft.  The  dominant  resonance  set  for 
a given  aircraft  are  distinct,  i.e.,  sets  for  other  air- 
craft, while  similar,  are  different. 

A prediction-correlation  processing  based  on  the  complex 
natural  resonances  of  a given  aircraft  can  be  used  to  iden- 
tify all  aircraft  whose  parameters  (dominant  complex  natural 
resonances)  are  known  a priori. 

Prediction-correlation  processing  of  echo  signal  data  for 
identification  of  aircraft  is  relatively  insensitive  to 
modest  amounts  of  noise  and  does  not  require  the  extraction 
of  complex  natural  resonances  from  an  unknown  (target)  tran- 
sient signal.  The  processing  therefore  is  applicable  to 
real  time  identification  of  fighter  aircraft.  In  this  fre- 
quency range  minor  substructure  changes  of  the  aircraft 
will  not  negate  the  identification. 
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Based  on  the  above  one  can  conclude  that  the  identification  of 
modern  fighter  aircraft  is  in-hand  provided  that  discrete  samples 
of  the  complex  echo  signal  spanning  the  frequency  range  2.0  to  24.0 
MHz  can  be  obtained.  Obviously  this  identification  method  is  not 
a retrofit  adjustment  of  existing  radar  systems.  It  is  stressed, 
however,  that  the  above  requirements  for  the  radar  are  specifically 
dictated  by  the  identification  task.  Stated  another  way,  given  com- 
plex echo  data  harmonically  spanning  the  2.0  to  24.0  MHz  frequency 
band  modern  fighter  aircraft  can  be  identified. 

A major  task  of  the  present  contracts  was  to  determine  if  modern 
fighter  aircraft  could  be  identified  using  similar  methods  if  the 
only  available  complex  echo  data  were  from  a significantly  higher 
frequency  band,  i.e.,  frequencies  above  24  MHz.  From  a complex  natural 
resonance  viewpoint  the  logical  such  substructure  resonances  are  those 
possibly  associated  with  the  stabilizer  resonances  in  the  30.0  MHz 
to  *56.0  MHz  range  and  cavity  resonances  nearer  the  upper  end  of  the 
band.  Neither  calculated  nor  measured  complex  echo  data  are  avail- 
able for  fighter  aircraft  in  the  30.0  MHz  to  TOO  MHz  frequency  region. 
Furthermore,  existing  experimental  facilities  which  measure  complex 
echo  data  at  harmonics  of  some  given  fundamental  frequency  are  not 
useful  because  the  frequency  increment  is  too  large.  To  obtain  the 
requisite  scattering  data  several  tasks  were  completed. 

7.  Preliminary  hand-operated  stepped  frequency  complex  echo 
measurements  were  made  in  the  30  MHz  to  50  MHz  region  using 
1/72  scale  models  of  the  F-104,  F-4  and  MIG-21  aircraft. 
Frequency  increment  for  these  data  was  5.0  MHz  and  a bi static 
angle  of  30°  and  vertical  (perpendicular  to  the  plane  of 
flight)  polarization  were  used.  The  measurements  were  made 
over  a range  of  aspect  angles  in  the  plane  of  flight.  Back- 
ground subtraction  at  each  discrete  frequency  was  used. 
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8.  From  the  measured  echo  data  contour  plots  of  the  magnitude 
of  the  echo  signal  for  frequencies  from  30.0  to  50.0  MHz 
and  aspect  angles  from  0°  to  180°  were  prepared.*  Both 

the  preliminary  measured  data  and  the  contour  maps  indicated 
that  30.0  MHz  to  56.0  MHz  was  the  frequency  region  where 
stabilizer  resonances  were  probable. 

9.  An  automated  stepped-frequency  measurement  range  was  in- 
strumented. The  facility  is  capable  of  rapidly  measuring 
and  calibrating  complex  echo  signals  at  stepped  frequencies 
from  (model  frequencies)  1.0  GHz  to  12.0  GHz.  The  system 
has  a variable  bistatic  angle  and  a variable  frequency  in- 
crement between  the  measured  stepped  frequencies.  Vector 
background  subtraction  is  used  at  each  measured  frequency. 
Accuracy  of  the  measurements  is  estimated  as  -1  dB  in  ampli- 
tude and  -10°  in  phase  at  each  frequency. 

’0.  Complex  echo  siqnal  data  in  the  frequency  range  30.0  MHz 
to  56.0  MHz  were  measured  for  six  aircraft  (F-104,  F-4, 

F15,  MIG-21,  F-104N,  F-105)  and  one  pencil-type  shape  with 
an  adjustable  cavity  length.  Certain  measurements  were 
also  made  in  the  frequency  range  from  56.0  MHz  to  ’06  MHz. 
Vertical  polarization  (perpendicular  to  the  plane  of  flight) 
were  used  for  transmitter  and  receiver. 

11.  Pulse  response  waveforms  with  sufficient  resolution  to  isolate 
substructure  features  of  the  aircraft  and  pencil  targets 
were  produced  from  the  stepped  frequency  data.** 


♦The  contour  maps  were  presented  in  a Quarterly  letter  dated 
13  September  1977. 

**The  measured  data  and  detailed  descriptions  of  the  targets  are  given 
In  Reference  24.  The  pulse  response  waveforms  are  also  given  In 
Reference  24. 


As  stated  in  the  footnote,  a much  more  detailed  discussion  of 
the  measurement  system,  the  measured  data,  the  simulated  pulse  response 
waveforms  and  the  targets  is  given  in  a data  report24.  We  concentrate 
here,  however,  on  certain  conclusions  most  pertinent  to  the  aircraft 
identification  problem. 

12.  In  the  frequency  range  from  30.0  MHz  to  56.0  MHz  for  ver- 
tical polarization  of  transmitter  and  receiver,  the  dominant 
contributor  to  the  echo  signal  from  modern  fighter  aircraft 
is  the  vertical  tail  stabilizer.  For  the  pencil  target 
at  frequencies  above  56.0  MHz  where  the  rectangular  cavity 
is  above  cutoff,  the  cavity  and  the  position  of  a short  with- 
in the  cavity  are  clearly  visable  in  the  echo  signal. 

’3.  Pulse  response  waveforms  from  aircraft  simulated  from  meas- 
ured complex  echo  signal  measurements  in  the  30.0  MHz  to 
56.0  MHz  range  did  not  demonstrate  the  stabilizer  assembly 
has  a high-0  complex  natural  resonance  in  this  frequency 
range.  A similar  conclusion  can  be  drawn  regarding  the 
cavity  at  frequencies  from  56.0  MHz  to  106.0  MHz  but  is 
more  tentative  because  the  measured  data  were  less  extensive. 

14.  In  most  cases,  the  dominance  of  the  vertical  stabilizer 
manifested  Itself  in  the  pulse  responses  as  a single  pulse 
in  the  response  waveform  for  a wide  variation  of  aspect 
angles. 

The  fact  that  the  stabilizer  assembly  on  fighter  aircraft  does 
not  have  a substructure  high-Q  complex  natural  resonance  in  the  fre- 
quency range  30.0  MHz  to  56.0  MHz  was  unexpected,  and  at  some  variance 
to  earlier  pulse  responses  obtained  from  echo  signals  at  frequencies 
below  ’A.O  MHz.  A completely  satisfactory  explanation  of  why  the 
tail  stabilizer  does  not  have  a high-Q  type  natural  resonance  in  this 
frequency  range  cannot  be  given  at  this  time. 


57 


15.  A method  for  combining  sampled  data  records  of  pulse  response 
waveforms  from  an  arbitrary  number  of  tarqet  aspects  has 
been  developed.  As  a result  It  Is  possible  to  find  that 
single  homogeneous  difference  equation  which  best  satisfies 
the  transient  records  of  a target  for  a variety  of  aspects. 

16.  It  has  been  demonstrated  that  an  elgenanalysis  solution 
for  the  coefficients  of  a linear  homogeneous  difference 
equation  is  more  versatile  than  Prony's  method  in  that  no 
a priori  assumptions  concerning  any  of  the  coefficients 
is  necessary. 

The  fact  that  only  a single  isolated  pulse  is  seen  in  the  tran- 
sient response  of  aircraft  as  simulated  from  complex  echo  data  in 
the  30.0  MHz  to  56.0  MHz  range  does  not  necessarily  preclude  use  of 
predictor-correlation  processing  for  identification  of  the  aircraft. 

A difference  equation  processing  can  be  applied  to  the  tail  (portions 
occuring  later  in  time  than  the  peak)  of  the  pulse  response  related 
to  the  stabilizer.  A similar  comment  can  be  made  concerning  the  tail 
of  the  oencil  target.  In  this  regard  one  can  just  as  easily  work 
with  the  envelope  of  the  pulse  response  since  it  is  not  anticipated 
that  the  complex  natural  resonances  deduced  from  the  coefficients 
have  any  physical  significance.  A physical  significance  for  the  coef- 
ficient is  a convenient  but  not  necessary  condition  for  Identification 
of  the  target. 

17.  It  has  been  demonstrated  that  a set  of  excitation  invariant 
difference  coefficients  can  be  found  for  the  pencil  target 
using  pulse  responses  from  various  aspects  as  simulated 
from  echo  signal  data  in  the  30.0  MHz  to  56  MHz  range.  It 
has  also  been  shown  that  a minor  structure  change  in  the 
vertical  stabilizer  yields  a different  set  of  excitation 
invariant  coefficients.  It  follows  that  aircraft  with  rela- 
tively similar  stabilizer  geometries  can  still  be  separately 


identified.  In  fact,  for  the  F-105  aircraft  model,  a set 
of  difference  equation  coefficients  with  stability  comparable 
to  the  pen  target  cases  has  been  found. 

Exploitation  of  the  developed  measurement  facility  and  processing 
of  data  obtained  from  this  facility  has  been  consistent  with  the  time 
and  funds  of  the  contract.  Fundamental  questions  concerning  the  use 
of  frequencies  above  24  MHz  for  the  identification  of  modern  fighter 
aircraft  have  been  answered.  Certain  apparent  misconceptions  con- 
cerning complex  natural  resonances  of  substructure  of  the  aircraft 
in  the  frequency  range  30.0  MHz  to  56.0  MHz  have  been  corrected.  Fi- 
nally a scheme  of  identification  of  fighter  aircraft  using  frequencies 
in  the  30.0  MHz  to  56.0  MHz  range  has  been  postulated  and  certain 
requisites  of  the  scheme  demonstrated. 

It  is  also  apparent  that  processing  tests  of  measured  data  already 
obtained  have  been  far  from  exhausted  and  that  many  new  measurements 
and  tests  have  been  suggested  by  the  present  results.  It  is  strongly 
recommended  that  research  programs  of  this  type,  which  explore  the 
echo  signals  of  modern  fighter  aircraft  at  frequencies  above  24  MHz, 
be  continued. 
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